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20 



30 



INPUT/OUTPUT 
UNIT 



DISPLAY 



40 



MAIN PROCESSING UNIT 

• IDENTIFY TARGET ELEMENTS 

• IDENTIFY SUBSET OF 
ACTIVE TARGET ELEMENTS AND 
CHOOSE POINT INSERT CRITERION 

• PREPROCESS BOUNDARY 
ACTIVE ELEMENTS 

• PROCESS REMIANING ACTIVE 
ELEMENTS PRODUCING A SET 
OF TERMINAL-EDGES 

• SELECT POINT OR POINTS FOR 
MESH INSERTION 

• RECALCULATE MESH 



MESH MEMORY 



INPUT GEOMETRY 




PROCESSING 
ELEMENTS 



TERMINAL-EDGES 



LONGEST-EDGE 
PARALLEL COMPUTER 
(FIG. 8) 



PROC. 1 



PROC. 2 



PROC. N 



50 



FIG.6 



MAIN 

PROCESSING 
UNIT(40) 



INPUT: - VERTICES OF OBJECT MODEL 

- BOUNDARIES OF OBJECT MODEL 

- INITIAL MESH 



200 



IDENTIFY TARGET-ELEMENTS ACCORDING AT LEAST 
ONE CHOSEN CRITERION BETWEEN: 

• BAD-SHAPED ELEMENTS 

• FINITE-ELEMENT ANALYSIS INDICATOR 

• USER-DEFINED SIZE FUNCTION 

• OTHER 
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- DEFINE SUBSET OF TARGET ELEMENT ST-EL 
(ONE OR MORE TARGET ELEMENTS) 
CHOOSE A POINT INSERTION CRITERION (DELAUNAY 
OR TERMINAL-EDGES-REFINEMENT), SET IFLAG 



220 




• PREPROCESS BOUNDARY ELEMENTS (FIG. 1 1) 

• RECALCULATE MESH (FIG. 10) 



PROCESS EACH ELEMENT T IN ST-EL 

• INITIALIZE EMPTY SUBMESH AND SET OF 
PROCESSING ELEMENTS WITH T 

• USE LONGEST-EDGE PARALLEL COMPUTER UNTIL 
SET OF PROCESSING-ELEMENTS IS EMPTY (FIG. 8) 

• SELECT ONE OR MORE POINTS ACCORDING 
CHOSEN POINT-INSERTION CRITERION (FIG. 9) 

• RECALCULATE MESH (FIG. 10) 



I 



240 



ACTUALIZE SUBSET ST-EL FROM TARGET -ELEMENTS 



270 



YES 




ACTUALIZE SET OF 
k TARGET ELEMENTS AND 
REPEAT FROM BOX 220 
OR EXIT 



FIG. 7 



LONGEST-EDGE PARALLEL COMPUTER (50) 



PROCESSOR i 
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PICK ELEMENT T OF LONGEST -EDGE LE 
FROM PROCESSING-ELEMENTS 
ADD T TO SUBMESH 



FIND SET OF ACTIVE-ELEMENTS 

• ELEMENTS NOT PRESENT IN SUBMESH 

• HAVING LE AS AN EDGE 

• HAVING LONGEST-EDGE > LE 



t 



FIND NON-ACTIVE ELEMENTS 

• ELEMENT NOT PRESENT IN SUBMESH 

• HAVING LE AS AN EDGE 

• HAVING LONGEST-EDGE <=LE 
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YES 



DATA IN 



NO 





| ACTIVE-FLFMENTS? 1 






450 






460 


ADD ACTIVE-F.I ,FMENTS 
TO PROCESSING-ELEMENTS 




ADD LE TO 
TERMINAL-EDGES 





r 



- ERASE ELEMENTS PRESENT IN SUBMESH 
470 FROM PROCESSING-ELEMENTS 

- EXIT 



PROCESSOR N 



FIG.8 



SELECT ONE OR MORE POINTS 

ACCORDING CHOSEN POINT INSERTION CRITERION 

(FROM BOX 240) 



YES 



510 



500 

IFLAG - DELAUNAY? 



NO 
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SELECT A POINT 
FOR DELAUNAY INSERTION 
(FIG. 12) 



SELECT POINTS FOR 
TERMINAL-EDGES REFINEMENT 
- MIDPOINTS OF TERMINAL-EDGES 




CONTINUE BOX 240 




FIG. 9 



RECALCULATE MESH 

(FROM BOX 235 OR 240) 



560 




POINT INSERTION 

ACCORDING DELAUNAY METHOD 

(EXTERNAL ROUTINE) 



TERMINAL-EDGES REFINEMENT 
FOR EACH TERMINAL-EDGE DO 
• BISECT ALL THE ELEMENTS 
SHARING TERMINAL-EDGE 
(INSERTION OF THE MIDPOINT 
OF TERMINAL-EDGE) 




OUTPUT: MESH 




FIG. 10 



PREPROCESS BOUNDARY ACTIVE ELEMENTS 

(FROM BOX 235) 



DO FOR EACH 
ELEMENT T IN 
SUBSET ST-EL 




• COMPUTE DB MAXIMUM LENGTH OF THE 
BOUNDARY EDGES OF T 

• COMPUTE DI MINIMUM LENGTH OF THE 
INTERIOR EDGES OF T 



640 




FIND POINT AUX 

EQUAL TO MIDPOINT OF THE 

LARGEST BOUNDARY EDGE OF T 



650 



• SELECT A BOUNDARY POINT 
IMPROVING THE BOUNDARY 
POINT DISTRIBUTION AROUND AUX 

• RECALCULATE MESH 

• ERASE T FROM SUBSET ST-EL 



t 

CONTINUE WITH NEX ELEMENT 
IN SUBSET ST-EL OR EXIT 



T REMAINS IN ~ n 
SUBSET ST-EL 



FIG. 11 



SELECT A POINT FOR DEL AUNAY INSERTION 

(FROM BOX 510) 



700 



- FIND G-EDGE THE GREATEST EDGE 
FROM TERMINAL-EDGES 

- FIND MPT MIDPOINT OF G-EDGE 



J30 

[S MPT A 
BOUNDARY POINT? 



740 



YES 



AUX-MPT 





NO 


COMPUTER DB DISTANCE FROM 
MPT TO THE BOUNDARY 







745 



750 

G-EDGE A BOUNDARY VERTI 
AND 
DB < K* G-EDGE? 



760 



YES 



AUX - PROJECTION 
POINT OF MPT 
OVER THE BOUNDARY 



NO 


. 770 


SELECT MPT 




(MPT IS INTERIOR VERTEX 


FAR FROM THE BOUNDARY) 



780 



SELECT A BOUNDARY POINT 

NEAR TO AUX IMPROVING 

THE BOUNDARY POINT DISTRIBUTION 




CONTINUE BOX 510 




FIG. 12 
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LONGEST-EDGE MESH DATA STRUCTURE 



EDGE: 






PV1,PV2 
LENGTH 


Pointers to the two vertices of EDGE 
Lenght of EDGE 




SI 
S2 


bet or pointers to tne elements naving 
longest- edge equal to EDGE 
Set of pointers to the elements having EDGE 
as an edge and with longest- edge greater than 
EDGE; S2 is empty for terminal- edges 


800 


FLFMENT: 






PV1.PV2, PV3, PV4 
PT1.PT2, PT3, PT4 

PLEDGE 


Pointers to the four vertices of ELEMENT 
Pointers to the four face-neighbor 
elements 

Pointer to the longest- edge of ELEMENT 
in the EDGE representation 


820 


VERTEX: 

X1,X2,X3 


Spatial coordinates of VERTEX 




VE-IND 


Vertex indicator value; equal to 0 if VERTEX 
belongs to the initial mesh; equal to the succesor 
of the maximum value between the VE-IND 
values of the vertices VI and V2, such that 
VERTEX was obtained by longest-edge bisection of 
the elements sharing edge VI -V2, otherwise 


840 


PGENEDGE 


Pointer to the generator edge of VERTEX, that 
is to preceding edge whose bisection produced 
the vertex VERTEX 




PELEMENT 


Pointer to one of the elements that share said 
VERTEX 




GEN-EDGE: 
PV1,PV2 


Pointers to the two vertices of the edge GEN-EDGE 
(this edge does not exist as an edge in the 
current mesh) 


860 


LENGTH 


Length of the edge GEN-EDGE 





FIG. 17 



INITIALIZATION OF THE LONGEST-EDGE MESH DATA STRUCTURE 



- FOR EACH VERTEX IN MESH DO 

CREATE VERTEX IN VERTEX REPRESENTATION; SET 
PGENEDGE - NULL 

- FOR EACH EDGE IN MESH DO 

CREATE EDGE IN EDGE REPRESENTATION; SET LENGTH, 
AND EMPTY SETS SI AND S2 

- FOR EACH ELEMENT IN MESH DO 

CREATE ELEMENT IN ELEMENT REPRESENTATION; SET 
POINTERS TO ITS 4 VERTICES AND ITS 4 FACE-NEIGHBOR 
ELEMENTS 



FOR EACH ELEMENT IN ELEMENT REPRESENTATION DO 
-FIND LONGEST -EDGE LE-EDGE AND SET PLEDGE (POINTER TO 
LE-EDGE IN EDGE REPRESENTATION) 
- FOR LE-EDGE IN EDGE REPRESENTATION DO 

• PUT T IN SET SI 

• PUT IN SET S2 ALL NEIGHBOR ELEMENTS SHARING EDGE 
LE-EDGE AND HAVING LONGEST-EDGE GREATER THAN 
LE-EDGE 



FOR EACH VERTEX IN VERTEX REPRESENTATION DO 

SET POINTER TO ONE OF THE ELEMENTS THAT SHARE VERTEX 
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FIG. 18 



MODIFICATION OF THE LONGEST-EDGE MESH DATA STRUCTURE 

BY VERTEX CREATION 

(FROM BOX 1150 IN FIG. 20) 

1000 



INPUT: - NEW VERTEX VN MIDPOINT OF EDGE OF VERTICES VI, V2 

- IFLAG (DELAUNAY OR TERMINAL-EDGE - REFINEMENT) 

- S-EL, SET OF ELEMENTS TO BE ELIMINATED 

- S-EDG, SET OF EDGES TO BE ELIMINATED 

- S NEW-EL, SET OF NEW ELEMENTS 

- SNEW-ED, SET OF NEW EDGES 



1020 

I 

-ELIMINATE ELEMENTS OF S-EL AND EDGES OF S-EDG FROM 
THE ELEMENT AND EDGE REPRESENTATIONS, RESPECTIVELY 

- CREATE VERTEX VN IN VERTEX REPRESENTATION, EDGES 
OF SNEW-ED IN EDGE REPRESENTATION, ELEMENTS OF 
SNEW-EL IN ELEMENT REPRESENTATION 

- ACTUALIZE INFORMATION OF OLD REMAINING EDGES AND 
VERTICES IN THE EDGE, ELEMENT AND VERTEX REPRESENTATION 



YES 




YES 

1060 



- CREATE GENEDE OF VERTICES VI, V2 IN GENEDGE REPRESENTATION 

- FOR VERTEX VN IN VERTEX REPRESENTATION DO 

- SET VE-IND- MAX (VE-IND VALUES OF VI AND V2) + 1 

- SET P GENEDGE EQUAL TO POINTER TO GENEDGE 



■4 

1080 



FOR VN, SET POINTER TO ONE OF THE ELEMENTS THAT 
SHARE VN 




FIG.19 



PARALLEL COMPUTER FOR TERMINAL-EDGE REFINEMENT 

(REPLACE BOX 570 IN FIG. 10) 



1100 

INPUT: MESH, SET OF TERMINAL-EDGES TO BE REFINED 
INITIALIZE ALL VERTICES AS NON-OCCUPIED AND 
TERMINAL-EDGES TO BE PROCESSED AS NON-PROCESSED 
FOR EACH TERMINAL-EDGE FIND S SET OF ALL ELEMENTS 
SHARING TERMINAL-EDGE, AND FIND SET SV OF NEIGHBOR 
VERTICES (VERTICES OF ELEMENTS OF S AND VERTICES OF 
ELEMENTS SHARING AT LEAST ONE VERTEX WITH ELEMENT 
OF S) 

CONTROL THE ASSIGNING OF TERMINAL EDGES TO 
INDIVIDUAL FREE PROCESSORS 

• SELECT TERMINAL-EDGE OF VERTICES VI, V2 AND 
MIDPOINT MPT 

• IF ALL VERTICES OF SV ARE NON-OCCUPIED MARK THEM 
AS OCCUPIED AND ASSIGN TERMINAL-EDGE, VI, V2, MPT, 
S AND SV TO FREE PROCESSOR 

• OTHERWISE SELECT NEW TERMINAL-EDGE 



I 



+ PROC. 1 



PROC. 2 



1150 



PROC. K 

- INPUT TERMINAL-EDGE OF VERTICES VI, V2 MIDPOINT MPT, 
AND ASSOCIATED SETS S AND SV 

• DIVIDE EACH ELEMENT OF S BY THE PLANE DEFINED BY MPT AND 
VERTICES OPPOSITE TO EDGE VI -V2 

• DEFINE SETS OF EDGES AND ELEMENTS TO BE ELIMINATED, AND 
SETS OF NEW EDGES AND NEW ELEMENTS TO BE CREATED 

• MODIFY MESH DATA STRUCTURE (FIG. 19) 

• MARK VERTICES OF SV AS NON-OCCUPIED 

• MARK TERMINAL-EDGE AS PROCESSED 

• FREE PROCESSOR 



PROC. N 



FIG. 20 



DEREFINEMENT METHOD 



INPUT: MESH 

VERTEX VX TO BE DEREFINED 



1200 



FIND NX, THE VE-IND VALUE OF NX 
GX, THE GENERATOR EDGE OF VX 



1220 



FIND VERTEX-SET ASSOCIATED WITHVX 
OF NODES (VZ, NZ, GZ) WITH 

VZ: NEIGHBOR VERTEX TO BE ELIMINATED 

NZ: THE VE-IND VALUE OF VZ 

GZ: THE GENERATOR EDGE OF VZ (FIG. 22) 



1240 



ORDER NODES OF VERTEX-SET 

IN DECREASING ORDER OF NZ AND, FOR EQUAL NZ, 

IN INCREASING ORDER OF GZ 



1260 



PARALLEL DEREFINEMENT OF THE ORDERED VERTICES 
OF VERTEX- SET 



1280 



YES 



1290 

EXISTS VX IN MESH? 




1300 



OUTPUT DEREFINED MESH 




FIG. 21 



VERTEX- SET CONSTRUCTION 

(FROM BOX 1040) 



1400 



INITIALIZE VERTEX- SET WITH NODE (VX, NX, GX) 
INITIALIZE VP-SET (SET OF NODES TO BE PROCESSED) 
WITH NODE (VX, NX, GX) 



1420 



PICK NODE (VY, NY, GY) FROM VP-SET 



1440 



FIND SET SN OF NODES (VZ, NZ, GZ), WHERE VZ IS 
NEIGHBOR VERTEX OF VY 



DO FOR EACH NODE 
(VZ, NZ, GZ) IN SN 




NO 



1480 




1520 



YES 



ADD NODE (VZ, NZ, GZ) TO VERTEX-SET 
ADD NODE (VZ, NZ, GZ) TO VP-SET 




FIG. 22 



PARALLEL COMPUTER FOR TERMINAL-EDGE DEREFINEMENT 

(FROM BOX 1280 IN FIG. 21) 



1630 



MESH 
MEMORY 
(BOX 10 
IN FIG. 6 
WITH 
LONGEST- 
EDGE 
DATA 

STRUCTURE 
OF FIG. 18) 



1600 

INPUT: SET OF VERTICES TO BE DEREFINED 

INITIALIZE ALL VERTICES AS NON-OCCUPIED AND VERTICES TO BE 
DEREFINED AS NON-PROCESSED 

CONTROL THE ASSIGNING OF ORDERED VERTICES TO INDIVIDUAL 
FREE PROCESSORS 

• SELECT VERTEX VA OF VE-IND VALUE NA AND GENERATOR EDGE 
VI -V2 

• FIND THE SET S OF NEIGHBOR VERTICES VZ OF VA SUCH THAT THE 
EDGE VA-VZ IS AN EDGE OF MESH 

• IF VERTICES VI, V2 BELONG TO S AND ALL VERTICES OF S ARE 
NON-OCCUPIED, MARK THEM AS OCCUPIED, MARK VA AS 
PROCESSIND AND ASSIGN VERTEX VA, VE-IND VALUE NA, 
GENERATOR-EDGE VI -V2 AND ASSOCIATED SET S, TO FREE 
PROCESSOR; OTHERWISE SELECT NEXT VERTEX 



PROC. 1 




PROC. 2 



1650 



PROC. K 

- INPUT: VERTEX VA, OF VE-IND VALUE EQUAL TO NA, GENERATOR 
EDGE OF VERTICES VI, V2 AND ASSOCIATED SET OF NEIGHBOR 
VERTICES S 

- FIND SET OF ELEMENTS TO BE ELIMINATED (SET OF ELEMENTS 
THAT SHARE VERTEX VA), AND SET OF EDGES TO BE ELIMINATED 

- FIND SET OF NEW ELEMENTS (ELEMENTS SHARING EDGE VI V2) AND 
SET OF NEW EDGES 

- MODIFY MESH DATA STRUCTURE (FIG. 24) 

- MARK VERTICES OF S AS NON OCCUPIED; MARK VA AS PROCESSED 

- FREE PROCESSOR 




PROC. N 



L 



FIG. 23 



MODIFICATION OF THE LONGEST-EDGE MESH DATA STRUCTURE 



BY DEREFIP 




1700 



INPUT: VERTEX VA TO BE ELIMINATED (MIDPOINT OF GENERATOR 
EDGE OF VERTICES VI, V2) 

- S-EL, SET OF ELEMENTS TO BE ELIMINATED 

- S-EDG, SET OF EDGES TO BE ELIMINATED 

- SNEW-EL SET OF NEW ELEMENTS 

- SNEW-ED, SET OF NEW EDGES 



- ELIMINATE ELEMENTS (OF S-EL) AND EDGES (OF S-EDG) FROM THE 
ELEMENT AND VERTEX REPRESENTATION, RESPECTIVELY 

- ELIMINATE VERTEX VA FROM VERTEX REPRESENTATION 

- CREATE EDGES OF SNEW-ED, AND ELEMENTS OF SNEW-EL IN EDGE 
AND ELEMENT REPRESENTATION, RESPECTIVELY 

- ACTUALIZE INFORMATION OF OLD REMAINING EDGES AND VERTICES 
IN THE EDGE, ELEMENT AND VERTEX REPRESENTATIONS 

- ELIMINATE GENARATOR EDGE OF VERTICES VI, V2 FROM GENEDGE 
REPRESENTATION 
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FIG. 24 



INTEGRATED MESH GENERATION METHOD 



1800 



INPUT: BOUNDARY MODEL 



i 



PRODUCE A SURFACE MESH OF THE OBJECT INCLUDING INTERFACES 
USE THE SURFACE -METHOD (IFLAG- DELAUNAY) TO PRODUCE A 
QUALITY SURFACE MESH ICLUDING INTERFACES 
CONSTRUCT A VOLUME DELAUNAY MESH USING SURFACE 
VERTICES 



i 



1810 



RESPECT GEOMETRY BOUNDARY AND INTERFACES 

- FIND SET S OF VOLUME ELEMENTS THAT INTERSECT BOUNDARY 
OR INTERFACES 

- FOR EACH ELEMENT T OF S AND UNTIL S IS EMPTY DO 

• FIND A SURFACE ELEMENT T-SURF OF THE SURFACE MESH THAT 
INTERSECTS T 

• USE THE SURFACE METHOD (WITH IFLAG - DELAUNAY) TO IMPROVE 
T-SURF (FIND A SURFACE POINT AND INSERT IT) 

• INSERT NEW POINT IN VOLUME-MESH 

• ACTUALIZE SET S (CURRENT VOLUME ELEMENTS THAT INTERSECT 
BOUNDARY OR INTERFACES) 



i 



1820 



IMPROVE THE VOLUME MESH USING IMPROVEMENT METHOD 
( IFLAG=DELAUNA Y) 



1830 



REFINE AND/OR DEREFINE MESH AND/OR MODIFY MESH GEOMETRY 
AND/OR VERTEX COORDINATES ACCORDING EXTERNAL 
(APPLICATION-RELATED) FUNCTIONS 
ACTUALIZE MESH MEMORY (FIG. 26) 



1840 



1870 



- CONSTRUCT DELAUNAY 
MESH OF CURRENT VERTICES 

-USE METHOD (IFLAG - 
DELAUNAY) TO 
IMPROVE MESH 




FIG. 25 



ACTUALIZE LOCALLY MESH DATA STRUCTURE 
ACCORDING VERTEX CHANGES 



(FROM BOX 1840) 



INPUT: SET S OF VERTICES CHANGED IN MESH AND THEIR 
COORDINATES 



1900 



FIND SET SV OF NEIGHBOR VERTICES CONNECTED AT LEAST 

WITH ONE VERTEX OF S 

SET S EQUAL TO THE UNION OF S AND SV 

FIND SET S-EDGE OF EDGES INCLUDING AT LEAST ONE VERTEX 
OF S 

FIND SET S-EL OF ELEMENTS INCLUDING AT LEAST ONE VERTEX 
OF S 



1920 



CONSTRUCT A RESTRICTED LONGEST-EDGE MESH DATA 
STRUCTURE USING THE PROCEDURE OF FIG. 18 WITH THE 
INITIAL MESH EQUAL TO THE SUBMESH (NOT NECESSARILY 
CONNECTED) DEFINED BY THE SETS S, S-EDGE AND S-EL 
ACTUALIZE THE GLOBAL MESH DATA STRUCTURE ACCORDING 
THE DATA OF THE RESTRICTED DATA STRUCTURE FOR 
SUBMESH, SET VE-IND=0 FOR THE VERTICES OF S 
ERASE GENERATOR EDGES HAVING VERTICES IN SV 
AND SET VE-IND - 0 FOR INVOLVES VERTICES 
ACTUALIZE VE-IND VALUES FOR THE NEIGHBOR VERTICES 
WHOSE VE-IND VALUES CHANGED TO 0 



1940 



1960 



OUTPUT: MODIFIED (NON-DELAUNAY) MESH 



FIG. 26 



